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Abstract 



We consider a family of Gagliardo-Nirenberg-Sobolev interpolation in- 
equalities which interpolate between Sobolev's inequality and the logarith- 
mic Sobolev inequality, with optimal constants. The difference of the two 
terms in the interpolation inequalities (written with optimal constant) mea- 
sures a distance to the manifold of the optimal functions. We give an ex- 
plicit estimate of the remainder term and establish an improved inequality, 
with explicit norms and fully detailed constants. Our approach is based on 
nonlinear evolution equations and improved entropy - entropy production 
estimates along the associated flow. Optimizing a relative entropy func- 
tional with respect to a scaling parameter, or handling properly second 
moment estimates, turns out to be the central technical issue. This is a 
new method in the theory of nonlinear evolution equations, which can be 
interpreted as the best fit of the solution in the asymptotic regime among 
all asymptotic profiles. 
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Consider the following sub-family of the Gagliardo-Nirenberg-Sobolev in- 
equalities 



with e = 9{p) := S-k 1 < p < 3^ if d > 3 and 1 < p < oo 



if d = 2. Such an inequality holds for any smooth function / with sufficient 
decay at infinity and, by density, for any function / G LP"'"-'^(IR'^) such that V/ 



1. Introduction and main results 



(1) 
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is square integrable. We shall assume that C^^ is the best possible constant 
in (1). In [16], it has been established that equality holds in (1) if f = Fp with 

(2) = (1 + VxGM"^ 

and that all extremal functions are equal to Fp up to a multiplication by a 
constant, a translation and a scaling. See Appendix A for an expression of C^^ . 
If d > 3, the limit case p = d/{d — 2) corresponds to Sobolev's inequality and 
one recovers the optimal functions found by T. Aubin and G. Talenti in [3, 23]. 
When p — 7- 1, the inequality becomes an equality, so that we may differentiate 
both sides with respect to p and recover the euclidean logarithmic Sobolev 
inequality in optimal scale invariant form (see [20, 25, 16] for details). 

It is rather straightforward to observe that Inequality (1) can be rewrit- 
ten, in a non-scale invariant form, as a non-homogeneous Gagliardo-Nirenberg- 
Sobolev inequality: for any / G L^^^ nV^'"^' 



(3) / \Vf\^dx+f \f\'+Ux>KpA( 

with 

The optimal constant Kp,i can easily be related with C^^ . Indeed, by op- 
timizing the left hand side of (3) written for f\{x) := X^/^'^p^ fi^x) for any 
X € W^, with respect to A > 0, one recovers that (3) and (1) are equivalent. 
The detailed relation between Kp^ and C^^ can be found in Section 7. 

Define now 

2(p-i) -n/d-p(d-4) 



Cm:=[j^) , M.:=y (l + |x|^) - = vr. 

Vp-i 

Consider next a generic, non-negative optimal function. 



1 



fMlj^)--=^"~'[CM + -\x-y\-'] VxG 



and let us define the manifold of the optimal functions as 

(p) 

We shall measure the distance to 9Hj^ with the functional 

7^(P)[/]:= inf / [g^~P {\f\^P - g^P) - ^^[\f\P+^ - gP+^)] dx . 
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To simplify our statement, we will introduce a normalization constraint and 
assume that / G L^^'(M^, (1 + \x\'^) dx) is such that 

fr,\ /Rri l/P^ dx _ d{p-l)a,M2~'' . x ._ ( , d+2~p(d~2) \ d-pl^d-A) 

^ ' {y dxf ~ d+2~p[d-2) ' ^* W - ^4 

Such a condition is not restrictive, as it is always possible to cover the general 
case by rescaling the inequality, but significantly simplifies the expressions. As 
we shall see in the proof, the only goal is to fix a = 1. 
Our main result goes as follows. 

Theorem 1. Let d > 2, p > 1 and assume that p < d/{d — 2) if d > 3. 
For any f £ LP+i nP^'2(R'^) such that Condition (5) holds, we have 

(7^(p)[/])' 



/ \Vf\'dx+f \f\P+Ux-Kp,d(f l/l'^rfx) >Cp,d 
where 7 is given by (4). 



iUfV^dxY 



The constant Cp^d is positive and explicit. We do not know its optimal value. 
See Appendix A for an expression of Cp^d, which is such that 

lim C„ d = and lim C„d = ^ ■ 

P^l+ p->d/{d~2)- 

The space 

LP+i p pi>2(]Rd) is the natural space for Gagliardo-Nirenberg in- 
equalities as it can be characterized as the completion of the space of smooth 
functions with compact support with respect to the norm || • || such that 



11^/112 + ll/llp+i- -'-'^ paper, we shall also use the notations 



15 

i„ - 1^1' l/l' -i^)"'. » "»t iiJ II, - «i no,,- 

Under Condition (5), we shall deduce from Theorem 4 that 

(6) 7^(^')[/]>CcK ll/lg^^"'^ inf_J/|2P-52p 



with 5 = d + 2— p{d + Q) for some constant Cck whose expression is given 
in Section 3, Eq. (13). Putting this estimate together with the result of Theo- 
rem 1, with 

^p,d '■= (^d,p CcK^ , 
we obtain the following estimate. 

Corollary 2. Under the same assumptions as in Theorem 1, we have 



/ |V/|2dx+ / \frUx-Kp,d(f l/l 

jRd. J^d \Jfid 



2p 



dx 



> \\f&-'' inf 
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The critical case p = d/{d—2) corresponding to Sobolev's inequality raises 
a number of difficulties which are not under control at this stage. However, 
results which have been obtained in such a critical case, by different methods, 
are the main motivation for the present paper. 

In [9, Question (c), p. 75], H. Brezis and E. Lieb asked the question of 

(p) 

what kind of distance to QJl^ is controlled by the difference of the two terms 
in the critical Sobolev inequality written with an optimal constant. Some 
partial answers have been provided over the years, of which we can list the 
following ones. First G. Bianchi and H. Egnell gave in [5] a result based on 
the concentration-compactness method, which determines a non-constructive 
estimate for a distance to the set of optimal functions. In [15], A. Cianchi, 
N. Fusco, F. Maggi and A. Pratelli established an improved inequality us- 
ing symmetrization methods. Also see [14] for an overview of various results 
based on such methods. Recently another type of improvement, which relates 
Sobolev's inequality to the Hardy-Littlewood-Sobolev inequalities, has been 
established in [17], based on the flow of a nonlinear diffusion equation, in the 
regime of extinction in finite time. Theorem 1 does not provide an answer 
in the critical case, but gives an improvement with fully explicit constants in 
the subcritical regime. Our method of proof enlightens a new aspect of the 
problem. Indeed, Theorem 1 shows that the difference of the two terms in 
the critical Sobolev inequality provides a better control under the additional 
information that ||/||2 2p finite. Such a condition disappears in the setting 
of Corollary 2. 

In this paper, our goal is to establish an improvement of Gagliardo- 
Nirenberg inequalities based on the flow of the fast diffusion equation in the 
regime of convergence towards Barenblatt self-similar profiles, with an explicit 
measure of the distance to the set of optimal functions. Our approach is based 
on a relative entropy functional. The method relies on a recent paper, [19], 
which is itself based on a long series of studies on intermediate asymptotics of 
the fast diffusion equation, and on the entropy - entropy production method 
introduced in [4, 2] in the linear case and later extended to nonlinear diffu- 
sions: see [21, 22, 16, 12, 11]. In this setting, having a finite second moment is 
crucial. Let us give some explanations. 

Consider the fast diffusion equation with exponent m given in terms of 
the exponent p of Theorem 1 by 

(7) p = <?=^ m = . 

^ ^ ^ 2m-l 2p 

More specifically, for m £ (0, 1), we shall consider the solutions of 

(8) ^ + V- [n (7?Vn™-i-2x)] =0 t>0, xGR'^ 
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with initial datum u{t = 0, •) = uq. Here is a positive parameter which 
does not depend on t. Let Uoo be the unique stationary solution such that 
M = J^d u dx = /jgd Uoo dx. It is given by 

Moo(a;) = (-?^ + - kp)""' yxeM."^ 

for some positive constant K which is uniquely determined by M. The follow- 
ing exponents are associated with the fast diffusion equation (8) and will be 
used all over this paper: 

d-2 d-1 , _ d 

rric := — ; — , mi := — ; — and mi 



d d d + 2 

To the critical exponent 2p = 2d/[d — 2) for Sobolev's inequality corresponds 
the critical exponent mi for the fast diffusion equation. For d > 3, the con- 
dition p £ {l,d/{d — 2)) in Theorem 1 is equivalent to m € {mi, 1) while for 
d = 2, p G (1, do) means m € (1/2, 1). 

It has been established in [21, 22] that the relative entropy (or free energy) 

J"[u|noo] := I [n™ - -u™ - m {u - Uoo)l dx 

m — I jRd L J 

decays according to 

^T[u{-,t)\Uoo] = -I[u{-,t)\Uoo] 

if w is a solution of (8), where 

m 



I[u{-,t)\uoo] := r] ■ 



2 

dx 



1 — m 

is the entropy production term or relative Fisher information. If m S [mi, 1), 
according to [16], these two functionals are related by a Gagliardo-Nirenberg 
interpolation inequality, namely 

(9) J"[u|tioo] < ^X[n|uoo] • 

We shall give a concise proof of this inequality in the next section (see Re- 
mark 1) based on the entropy - entropy production method, which amounts to 
relate ^I[u{-,t)\uoo] and X[ii(-, t)|uoo]- We shall later replace the diffusion pa- 
rameter r] in (8) by a time-dependent coefficient cr(t), which is itself computed 
using the second moment of u, J^a \x\'^ u{x, t) dx. By doing so, we will be able 
to capture the best matching Barenblatt solution and get improved decay rates 
in the entropy - entropy production inequality. Elementary estimates allow to 
rephrase these improved rates into improved functional inequalities for / such 
that l/P^ = u, for any p € {l,d/{d — 2)), as in Theorem 1. 

This paper is organized as follows. In Section 2, we apply the entropy - 
entropy production method to the fast diffusion equation as in [11]. The key 
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computation, without justifications for the integrations by parts, is reproduced 
here since we need it later in Section 6, in the case of a time-dependent diffusion 
coefficient. Next, in Section 3, we estabhsh a new estimate of Csiszar-Kullback 
type. By requiring a condition on the second moment, we are able to produce 
a new estimate which was not known before, namely to directly control the 
difference of the solution with a Barenblatt solution in L^(M'^). 

Second moment estimates are the key of a recent paper and we shall 
primarily refer to [19] in which the asymptotic behaviour of the solutions of 
the fast diffusion equation was studied. In Section 4 we recall the main results 
that were proved in [19], and that are also needed in the present paper. 

With these preliminaries in hand, an improved entropy - entropy produc- 
tion inequality is established in Section 5, which is at the core of our paper. 
It is known since [16] that entropy - entropy production inequalities amount 
to optimal Gagliardo-Nirenberg-Sobolev inequalities. Such a rephrasing of our 
result in a more standard form of functional inequalities is done in Section 6, 
which contains the proof of Theorem 1. Further observations have been col- 
lected in Section 7. One of the striking results of our approach is that all 
constants can be explicitly computed. This is somewhat technical although 
not really difficult. To make the reading easier, explicit computations have 
been collected in Appendix A. 



2. The entropy - entropy production method 

Consider a solution u = u{x, t) of Eq. (8) and define 
z{x,t) := r]Vu"^~^ -2x 
so that Eq. (8) can be rewritten as 

^ + v.(..) = o. 

To keep notations compact, we shall use the following conventions. If ^ = 
{Aij)f j^i and B = {Bij)f j^^ are two matrices, lai A : B = J2tj=i "^ij ^ij ^""-^ 
I^P = A : A. If a and b take values in M.'^, we adopt the definitions: 

a-b = 2_^aibi, V-a = 2^— , a(g)6 = (oj 6^)^ ,,=1 , V(^a=\j-^j 

Later we will need a version of the entropy - entropy production method 
in case of a time-dependent diffusion coefficient. Before doing so, let us recall 
the key computation of the standard method. With the above notations, it is 
straightforward to check that 

^ =ry(l -m)V(u"'~2V- (nz)) and V (g) z = ?? V Vn""^ - 2Id . 
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With these definitions, the time-derivative of r/X[ii|noo] = /]gdti|-zp dx 
can be computed as 

d f I ,2 , f du , ,2 , ^ f dz 
— u\z\ ax = — \z\ ax + 2 uz ■ — ax . 

at JiRd Jffid ot JiRd ot 

The first term can be evaluated by 



ot 



V ■ {uz) \z\'^dx 
= 2/ uz® z :V (i$> z dx 

= 2r] uz ® z -.V (g)Vu'^~^ dx - 4: u\z\'^dx 
= 2T]{l-m) u"'~'^Vu(E)V ■■ {uz (g) z) dx - A u\z\'^dx 
= 27]{l-m)[ u"'~^{Vu- zfdx+ 2r]{l-m) [ n™"^ (Vu • z) (V • z) 
+ 2i]{l-m) u™"^ (z (g) Vn) : (V (g) z) - 4 / u\z\'^dx . 

The second term can be evaluated by 

f , 
2 uz ■ — dx 
dt 



= 2r]{l-m) (uz • V) V • (liz)) dx 

JRd ^ ' 

= -2r]{l-m)f u""-^ (V ■ {uz)Ydx 

= -27j{l-m)[ \u'^{V ■ zf + 2u'^-\Vu ■ z) {V ■ z) + u'^'-^Vu ■ zf] dx . 

jRd 

Summarizing, we have found that 
[ — |zpdx-|-4 / tijzpdx 

JRd ot JRd 

= -2r] (1-m) / \u^{V-zf+u {Vu ■ z) {V ■ z)-u {z(S)Vu) : (Vfg) z)] dx . 

Using the fact that 

d'^z^ d'^z^ 



dxi dxj dx"^ ' 
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we obtain that 

/ u"^-^ [Vu • z) (V • z) dx 



m jRd m jRd oxi ox 



and 



^'""^(z^Vu) :{V(g)z) dx 

1/" ^i„i9. If ^ ^ — > 4 d /. 

= ~ / 

m jRd 



I vJ^\Vz\'^ dx + — I u'^Yz'^dx 
JRd m jRd /-^ dx^ 



i,j=l -^-^3 



can be combined to give 

■"""^ [ u iVu ■ z) (V • z) - n Vu (g) z : V (Xi z] dx 



u 



-— I -u™ (V -zf dx + — f n*" |Vz|2 dx 
m Jud m jRd 



This shows that 



d 



(10) ^/ u\z\'^dx + 4: u\z\'^ dx 

dt jRd J^d 

1 — m 



2rj -f u^'ilVz]^ -{1-m) {V -zf) dx . 

m Jud ^ ' 



By the arithmetic geometric inequahty, we know that 

|Vz|2 - (1 - m) (V • zf > 

if 1 — m < that is, if m > mi. Altogether, we have formally established 
the following result. 

Proposition 3. Let d > 1, m £ {mi, 1) and assume that u is a non- 
negative solution of (8) with initial datum uq in L^{R'^) such that and 
X I—)- l^ptiQ are both integrable on W^. With the above defined notations, we 
get that 

^I[u{-,t)\u^] < -AI[u{-,t)\u^] yt>o. 

The proof of such a result requires to justify that all integrations by parts 
make sense. We refer to [12, 13] for a proof in the porous medium case (m > 1) 
and to [11] for mi < m < 1. The case m = 1 was covered long ago in [4]. 

Remark 1. Proposition 3 provides a proof of (9). Indeed, with a Gronwall 
estimate, we first get that 

I[u{-,t)\uoo]<I[uo\uoo]e-^* yt>0 
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ifl[uo\uoo] is finite. Since I[u(- , t)\uao] is non-negative, we know that 

lim I[u{-,t)\uoo] = , 

which proves the convergence of u{-,t) to u^q as t ^ oo. As a consequence, we 
also have limt^oo J^[u{- , t)\uoo] = and since 

^{l[u{-,t)\Uoo] - 4:T[u{-,t)\Uoo]) = ^I[u{-,t)\Uoo] + AI[u{-,t)\Uoc] < , 

an integration with respect to t on (0, oo) shows that 

I[uo\Uod] - iT[uo\Uoo] > , 
which is precisely (9) written for u = uq. 

3. A Csiszar-Kullback inequality 

Let m G {rhi, 1) with mi = and consider tiie relative entropy 

[u] := / [n™ -B^-m B^~' (n - B^)] dx 

771-1 jRd L J 

for some Barenblatt function 

(11) B4x)■.= a~^CM + ^\xf)^' VxeM"^ 

where cr is a positive constant and Cm is chosen such that = M > 0. 

With p and m related by (7), the definition of Cj\/ coincides with the one of 
Section 1. See details in Appendix A. 

Theorem 4. Let d > 1, m £ {rhi, 1) and assume that u is a non-negative 
function in L^(R'^) such that n"* and x ^ \x\'^ u are both integrable on M'^. // 
= M and Jj^d \x\'^ u dx = J^d |xp i?o- dx, then 

4^ > ^rZm, \\u-B^\, + ^ I \X\^\U- B^\ dx)' . 

Notice that the condition j^a \x\'^ u dx = j^d \x\'^ B(j dx is explicit and 
determines a uniquely: 

a = — — I \x\^ u dx with Km '■= [ Ixl"^ Bi dx . 

For further details, see Lemma 5 and (20) below, and Appendix A for detailed 
expressions of Km and J^d -Bf* dx. With this choice of a, since B^~^ = 
Cm + cjf we remark that J^d B'^~^ (n - B^) dx = so 

that the relative entropy reduces to 

T,[u]:=^—- f [u"'-B^]dx 
m - 1 jRd 
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Proof of Theorem 4. Let v := u/B^ and dii„ := dx. With these nota- 
tions, we observe that 

/ {v-l) dfi^ = f (n - B^) dx 

Cm f {u - B^) dx + a^^""^-"^^ f jxp (u - 5^) = . 



= era 
Thus 



(l-m) 



{v - 1) dfia = {v - I) dfi^ - (1 - dl^la = , 

' Jv>l Jv<l 

which, coupled with 

/ {v-l)d^ia+ {'^-v)d^„=j \v-l\diia, 

Jv>l Jv<l JR-i 

imphes 

/ \u- B^\B'l'-^ dx = ( \v-l\dfi„ = 2[ \v - 1\ dfi^ . 
Jw^ Jm.'i Jv<i 

On the other hand, a Taylor expansion shows that 

^Au] = / \v^-l-7n{v- 1)1 dfi^='^ [ \v - 1|2 dpi, 

for some function ^ taking values in the interval (min{l,u},max{l,t)}), thus 
giving the lower bound 

I J v<\ ^ Jv<l 

Using the Cauchy-Schwarz inequality, we get 

= / \v -l\Ba' Ba' dx] < \v-l\'^dn, B^dx 

and finally obtain that 

^ r , m (/^<i \v - 1| di^L^Y _ m (/igrf \u - B,\ B]^-^ dxf 
"^''^-2 yB-^dx 8 yB^dx 

which concludes the proof. □ 

Notice that the inequality of Theorem 4 can be rewritten in terms of 
|/|2p = u and = B, with p = 1/(2 m — 1). See Appendix A for the compu- 
tation of /jgd dx, cj. Cm and Km in terms of f^d |xp u dx and M^,. In the 
framework of Corollary 2, we observe that Condition (5) can be rephrased as 

1 ll^l|2p 1 Il/ll2,2p 

Altogether we find in such a case that 
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with 
(13) 



Cr-t^- - El± d+2-p(d~2) J-hw A,,rl-7 



p+1 32 p 

Remark 2. Various other estimates can be derived, based on second order 
Taylor expansions. For instance, as in [16], we can write that 



with V := u/B(j and ip{s) := s^^^, and get 



m 



-2m 



1 



Ll/">(Rd,rf^^) 



Using |b'"||Li/™(Rd,dM.) 
In'" - dx = 



|l|lLV'"(IRd,(i/^^)} 
II Dmlim ; 

< IIt;™ - 1 

62/ t/ie Cauchy-Schwarz inequality, we find 



II pm||l-m 
d^J.^) 11-° o- 111 



B. 



cr 111 



m 2 



2 m 



IS. 



o- 111 



With f = 2 J t/iis a/so ^ii^es another estimate of Csiszdr-Kullback type, 
namely 

Kp^d ,,,„ii „ii2 



f (p-i) 

2,2p 



, inf 

i(d+2-p(d-2)) (p) 

2p y d 



I/I 
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p+i 



/or some positive constant Hp^d^ which is valid for any p € (l,oo) if d = 2 
and any p G (1, -^j^] if d > 3. Also see [24, 12, 10, 18] for further results on 
Csiszdr-Kullback type inequalities corresponding to entropies associated with 
porous media and fast diffusion equations. 

4. Recent results on the optimal matching by Barenblatt solutions 

Consider on M'^ the fast diffusion equation with harmonic confining po- 
tential given by 
du 



(14) 



dt 



+ V 



u 



X 



t > , x£ 



with initial datum uq. Here cr is a function of t. Let us summarize some results 
obtained in [19] and the strategy of their proofs. 

Result 1. At any time t > 0, we can choose the best matching Barenblatt as 
follows. Consider a given function u and optimize A 1— )• -Fa [it]- 
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Lemma 5. For any given u S LY(M'^) such that u'" and are both 

integrahle, if m £ (mi,l), there is a unique A = A* > which minimizes 
A I— 7- J^A [u] , and it is explicitly given by 

A* = — — / Ixl"^ u dx 



Km 

where Km = |xp Bi dx. For A = A*, the Barenblatt profile B\ satisfies 

\x\'^ B\ dx = I \x\'^ u dx . 



As a consequence, we know that 

Of course, if u is a solution of (14), the value of A in Lemma 5 may depend 
on t. Now we choose a{t) = A(t), i.e., 

(15) a{t) = -^[ \x\'^u{x,t) dx yt>0. 

Km jr'^ 

This makes (14) a non-local equation. 

Result 2. With the above choice, if we consider a solution of (14) and compute 
the time derivative of the relative entropy, we find that 

However, as a consequence of the choice (15) and of Lemma 5, we know that 



^TM) =0, 



da 



and finally obtain 



dx . 



(16) - J-,(,)[n(.,0] = Y / ^ ^ 

dt ^ ' 1 — m 

The computation then goes as in [7, 8] (also see [21, 22, 16] for details). With 
our choice of o", we gain an additional orthogonality condition which is useful 
for improving the rates of convergence (see [19, Theorem 1]) in the asymptotic 
regime t — )■ oo, compared to the results of [8] (also see below). 

Result 3. Now let us state one more result of [19] which is of interest for the 
present paper. 

Lemma 6. With the above notations, if u and a are defined respectively 
by (14) and (15), then the function t i— ?• a{t) is positive, decreasing, with 
(Too := limt^oo cr{t) > and 

(17) a' it) = -2 d iljipl ("c) ^ [^(.^ < . 

m j\m 
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The main difficulty is to establish that cJoo is positive. This can be done 
with an appropriate change of variables which reduces (14) to the case where a 
does not depend on t. In [19], a proof has been given, based on asymptotic 
results for the fast diffusion equation that were established in [16, 7, 6, 8]. An 
alternative proof will be given in Remark 3, below. 

5. The scaled entropy - entropy production inequality 

Consider the relative Fisher information 



{m—Tric) 



m 



dx . 



1 — m . 

By applying (9) with Uoo = Bi and = 1 to x i— )■ cr'^/^ u{y/a x) and using the 
fact that Bi{x) = a^^'^ B„{^/a x), we get the inequality 

Now, if a is time-dependent as in Section 4, we have the following relations. 

Lemma 7. If u is a solution of (14) with a{t) = J^d \x\'^u{x,t) dx, 
then a satisfies (17). Moreover, for any t > 0, we have 



(18) 
and 



dt 



(19) Ax,(,)[^(.,t)]< 



4 H — [m — nic) (m — mi) a — -— 
2 a[t) 



Proof. Eq. (17) and (18) have already been stated respectively in Lemma 6 
and in (16). They are recalled here only for the convenience of the reader. It 
remains to prove (19). 

For any given a = (j{t), Proposition 3 gives 

|x,(,)[<,t)] = (^I,[u{;t)]) +a'{t)(^I,[u]) 

dt \dt ^|A=<7{t) ^dX ^|A=(7(t) 



<-4I„^,^[u{;t)]+a'{t)( — Ix[u] 



d 



\X=ait) 



Owing to the definition of Ix, we obtain 

^2:a[u] = ^ (m - nic) jlx[u] 

_ Af / 2u (Vn"^-i - VS5"-i) • — (VB 

1-m JRi ^ ^ ^ dX^ 

By definition (11), VS^'^x) = 2xA-i("'-"=\ which implies 



dx 



dX ^ ^ 



- — (m — mr) X . 
X ^ ' 
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Substituting this expression into the above computation and integrating by 
parts, we conclude with the equahty 



^2:aM = ^ (m - rUc) jlx[u] 



2d 

— (m - mc) 



2 7nA^f('"-™=) 



\x\ u dx — d ti™ dx 



1 — m 

A simple computation shows that 

(20) df B''^dx = -[ x-VBY'dx = ^^f \x\^ dx = Km 

jRd J^d i-m Jm.<i l-m 

and, as a consequence, if A = cr = J^d jxp u dx, then 

/ \x\'udx = d BTdx, 

l-m JM.'i 

and finally 

(21) ^^^N " ^ " - 4d(l - m)J^x[u]) . 
Altogether, we have found that 

^(X,(t)[tx(-,t)]) + 4I^^t)[u{;t)] 



<^{m- mc) ^-j^ [n] - 4 d (1 - m) T^^t) M 



^ 2 '"^ cj(t) 

The last term of the right hand side is non-positive because by (17) we know 
that cr'(t) < and 

^a(t) M - 4 d (1 - m) [u] 

= d{l-m)(l^(^t)[u\- 4-7v(t)[n]) + d (m - mi)X^(t) [u] 

> d{m- mi)Z„(t)[u] > . 

This implies (19). □ 

To avoid carrying heavy notations, let us write 

fit) := T,^t)H;t)] and j{t) := J^(t)H;t)] 

and denote /(O), j(0) and a{0) respectively by /o, jo and fio. Estimates (17), 
(18) and (19) can be rewritten as 

/' = -i < 

(22) { ct' = -Kifji^™-™^)/ < 
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where the constants Ki, i = 1, 2, are given by 

{\ _ jYiX^ \ 

Ki := 2d — and K2 := — (m — iJic) (m — mi) S . 

mKu 2 

Using the fact that hmf_j.oo /(i) = hmt_j.oo /(i) = 0, as in the proof of Propo- 
sition 3, we find that i{t) - Af{t) > and f{t) < fo e"^* for any t > 0. 

Remark 3. The decay of a can be estimated by 

(af (1-)) = ^ (1 - m) / < ^ (1 - m) fo e"^* , 

thus showing that a^J^^ '"^ > ctq^^ ™^ — ^47^^~ /o- Since uq and B^^ have 
the same mass and second moment, we know that 

fo = z / [B -UQ)dx= — "2 (T§ - r 

1 — m J^d ^ " ^ a(l— m)^ (1 — mj 

Hence we end up with the positive lower bound 

Z 4 m A A/ JRd 

Prom (22) we get the estimates a{t) < ctq for any t > and 

j' -Af'=j'+Aj<K2j- = Kl K2 (1—) ff'<Kl K2 Ct"^ ('"'"^ / /' 

a 

Integrating from to oo with respect to t and taking into account the fact that 
hmt^oo fit) = hmt->oo f{t) = 0, we get 

A I- ^ 1 -f (1-™) ^2 

- Jo + 4 /o < - - Ki K2 2 . 

By rewriting this estimate in terms of [uo] = fo, TctqIuq] = jo and after 
omitting the index 0, we have achieved our key estimate, which can be written 
using 

Cm,d ■■= 7, 7^ ("i - "^c) {m - mi) (1 - m)^ 

2m Km 

as follows. 

Theorem 8. Let d > 1, m £ (mi, 1) and assume that u is a non-negative 
function in L^(M'^) such that n™ and x i— ?• |xp u are both integrable on M'^. Let 
a = /jgd |xp u{x) dx where M = J^a u{x) dx. Then the following inequality 
holds 

(23) ATM + C^4 ^^^1^ < Uu] ■ 

Recall that Km = KiM^, with 7 = -^^^^i^-y- See Appendix A for details. 
Notice that this definition of 7 is compatible with the one of Theorem 1 if 
p = l/(2m - 1). 
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Remark 4. If we do not drop any term in the proof of Proposition 3 and 
Lemma 7, an ODE can be obtained for j , based on (10) and (21) and we can 
replace (22) by a system of coupled ODEs that reads 

f = -j<o 

a' = -Ki / < 

/ + ^3=d{l-m){j-Af)^ + K2j^ -r 

where r := 2f7'^(™-™-) /jg^ u'" (|Vz|2 - (1 - m) {V ■ zf) dx > and z := 

It is then clear that the estimates o" < co and j' + 4j < K2j which 
have been used for the proof of Theorem 8, are not optimal. 



6. Proofs of Theorem 1 and Corollary 2 

Let us start by rephrasing Theorem 8 in terms of / = u™~^/^. Assume 



that 



M=f udx= I \f\^P dx and a = f \x\^ u dx = f \x\'^ \f\'^P dx 

where p = 1/(2 m — 1) and using the notation fj^\ ^ £ 571 J'^ defined in Sec- 
tion 1, consider the functional 



dx 



In preparation for the proof of Theorem 1, we can state the foUowing resuh. 

Corollary 9. Let d > 2, p > 1 and assume that p < d/{d — 2) if d>3. 
For any f G LP+i nPi'2(M'^) such that Condition (5) holds, we have 



f \Vf\'dx+f \f\P+Ux-KpA[ |/|2Pd^)^>Cp,,-i^^^^ 



where 7 is given by (4). 



This results is slightly more precise than the one given in Theorem 1, as 

(p) 

we simply measure the distance to a special function in , the one with 
same mass and second moment, centered at 0. The constant Cp^d is the same 
as in Theorem 1: see Appendix A for its expression. 
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Proof. By expanding the square in I„ [u] and collecting the terms with the 
ones of J> [u] , we find that 

+ f y-- dx + j^imKuai^'-^^ - [ B^dx). 

The last term of the right hand side can be rewritten as 

l-m cr J l~m {d+2) m-d 

with 7 = ^^^^'^ !^ (as in the previous Section) and Ci = mI (see Appen- 
dix A for details). Consequently Inequality (23) can be equivalently rewrit- 
ten as 

(24) cji^"'-'"^) J^^ I Vu™-^ |2 dx + d ^ J^^ dx 

— l—m {d+2jm—d ^ 

. d^ (m—mc) {m—mi) {1—m)^ t[u\) 

This inequality is invariant under scaling and homogeneous. As already no- 
ticed in (12), Condition (5) means a = a^, that is cjf = d 
Using the explicit expressions that can be found in Appendix A and reexpress- 
ing all quantities in terms of p = 2m-i completes the proof of Corollary 9. See 
Appendix A for an expression of Qp^d- ^ 

Proof of Theorem 1. It is itself a simple consequences of Corollary 9. 

Let us consider the relative entropy with respect to a general Barenblatt 
function, not even necessarilynormalized with respect to its mass. For a given 
function u € L3,_(M'^) with n"^ € L^(]R'^) and \x\'^ u G L^(R'^), we can consider 
on (0, oo) X X (0, oo) the function h defined by 

h{C, y, a) = [u"^ - i?^,,, - m 5""^ {u - Bc,y,.)\ dx 

where Bc,y,a is a general Barenblatt function 

Bc,yA^) := fJ-i (C + i |x - ^ ^xGR'^. 
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An elementary computation shows that 
dh maf^^-"') 



dC 1 — m 

-| (m-inc) 



u - Bc,y,a) dx , 

Vyh = / {x -y){u- Bc,y,a) dx , 

C (n - Bc,y,a) dx 

JR'' 



dh d N 
— =m-a 2 1™ 
da 2 



m — rrir 1 



m a 



/ \x - y\'^ {u - Bc,y,a) dx 



Optimizing with respect to C fixes C = Cm, with M = J^a u dx. Once 
C = Cm is assumed, optimizing with respect to a amounts to choose it such 
that Jj^d |xp Bc,y,a dx = /jgd \x — y\'^ u dx as it has been shown in Lemma 5. 

This completes the proof of Theorem 1, since 7^(p)[/] > R(P)[/] by def- 
inition of TZ^P^ (see Section 1). Notice that optimizing on y amounts to fix 
the center of mass of the Barenblatt function to be the same as the one of u. 
This is however required neither in the proof of Corollary 9 nor in the one of 
Theorem 1. □ 

Proof of Corollary 2. It is a straightforward consequence of Theorem 1 
and of the Csiszar-Kullback inequality (6) when / G is such that 

II/II22P is finite. However, ||/||2 2p does not enter in the inequality. Since 
smooth functions with compact support (for which ||/||2 2p is obviously finite) 
are dense P^'^(]R'^), the inequality therefore holds without restriction, by den- 
sity. □ 

7. Concluding remarks 

Let us conclude this paper with a few remarks. First of all, notice that 
Theorem 4 gives a stronger information than Theorem 1, as not only the 
L^{M.'^,dx) norm is controlled, but also a stronger norm involving the second 
moment, properly scaled. 

No condition is imposed on the location of the center of mass, which simply 
has to satisfy (J^dXudx)'^ < J^dU dx J^^d dx = aM Km according to 

the Cauchy-Schwarz inequality. Hence in the definition of TZ[f] and TZ^^[f] (in 
Theorem 1) as well as in Corollary 2, the result holds without optimizing on 
y G W^. In [8, 19], improved asymptotic rates were obtained by fixing the center 
of mass in order to kill the linear mode associated to the translation invariance 
of the Barenblatt functions. Here this is not required since, as i — >• 00, the 
squared relative entropy is simply of higher order. Our improvement is better 
when the relative entropy is large, and is clearly not optimal for large values of t. 
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Our approach differs from the one of G. Bianchi and H. Egnell m [5] 
and the one of A. Cianchi, N. Fusco, F. Maggi and A. Pratelh, [15]. It gives 
fuhy explicit constants in the subcritical regime. The norms involved in the 
corrective term are not of the same nature. 

Let us list a series of remarks which help for the understanding of our 
results. 

(i) Scaling properties of the Barenhlatt profiles. Consider the scaling X ^ u\ 
with u\{x) := A'^u(Ax) for any x G M"^. Then we have 

W I |2 W 1 W I |2 ^ ^ 

cja := -rr- / F u\ dx = -— / \x\ u dx = 
and may observe that 

B„,{x) = X'^B^iXx) . 
As a consequence, we find that = A''^'""^) 

(ii) Homogeneity properties of the Barenhlatt profiles. Similarly notice that for 

2 (1-m) 2(l-m) 

any m G (mi,l), we have Cm = CiM di^-mc) and Km = Ki M d(m-mc) , 
Let u\ := Xu and denote by Ba-^ the corresponding best matching Barenhlatt 
function. Using the fact that ||?iA|li = AM if \\u\\^ = M and observing that 

^ 2(i-m) r r 

KxM = Km X d{n,-mc) and / \x\'^ ux dx = X \x\'^ u dx , 

we find 

2 (1-m) 



ax = / ^^A dx = A'*('"-'"c) a . 

KxM 



2 (1-m) 



Since Cxm = A ^cm-mc) Cm, we find that 



d 



1 



/ 2 (1-m) I 2 (1-m) , ,2 \m-l 

B,^{x) = (^A^(— '»^) a\ I A ^(^-^c) Cm + 2I1L) = \B,{x) . 

\ xd(.m-mc) 

As a consequence, we find that Jv[ua] = X^ Fa[u]. 

(iii) The m = 1 limit. As m — )• 1, which also corresponds to p — )■ 1, we 
observe that the constant Cp^d in Theorem 1 has a finite limit. Hence we get 
no improvement by dividing the improved Gagliardo-Nirenberg inequality by 
{p — 1) and passing to the limit p — )• 1+, since TZ^^^[f] = 0{p — 1). By doing 
so, we simply recover the logarithmic Sobolev inequality as in [16]. 

This is consistent with the fact that, as m — )■ 1_, we have Cm,d ~ (l"""^)^; 
a = 0{K^j) = 0(1 — m) and, since 

dM \^ f d M , 



^ ^ 27: ]j;]2 ti J ^\ 2 f^d\x\'^ u dx^ 
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we also get that J>[n] ~ J^dU log (-^j dx. Hence, in Theorem 8, the addi- 
tional term in (23) is of the order of 1 — m and disappears when passing to the 
limit m — ?► 1_. 



Appendix A. Computation of the constants 

Let us recall first some useful formulae. The surface of the d — 1 dimen- 
sional unit sphere S is given by IS'^-i] = 2TT'^/'^/T{d/2). Using the integral 
representation of Euler's Beta function (see [1, 6.2.1 p. 258]), we have 



1 -|- X ax = 7r2 



Via) ■ 

With this formula in hand, various quantities associated with Barenblatt func- 
tions can be computed. Applied to the function B{x) := (l -|- |xp) '"^^ , x G M"', 



we find that 



-p/ d(m-me) \ 

M.:= / Bdx = 7rl-^^^^^ 



, l—m / 

Notice that when M = M*, B = Bi with the notation (11) of Section 3. As a 

consequence, for Bi{x) = {Cm + """^ j a simple change of variables shows 
that 

M := / Bidx= [Cm + \x\^)~ dx = C, 



which determines the value of Cm, namely 



2(l-"i) 



2(l-m) 2(l-m) 

A useful equivalent formula is Cm = CiM d{m-mc) where Ci = M*'' . 
By recalling (20) and observing that 

/ B"^ dx= [ B"^-^ Bidx= I {Cm + kp) Bi dx = M Cm + Km 

where Km ■= J^d\x\'^ Bi dx, using MCm = Ci with 7 = -^^7^1^^, we 

find that 

(25) 

Km= , CiM^ and / = -C^M^ . 

{d + 2)m-d jRd ^ ^d + 2)m-d 

Consider the sub-family of Gagliardo-Nirenberg-Sobolev inequalities (1). 
It has been established in [16, Theorem 1] that optimal functions are all given 
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by (2), up to multiplications by a constant, translations and scalings. This 
allows to compute C^^ . All computations done, we find 

^GN _ f {P~ir+^ Y ( d+2-p{d-2) \T? ( ^^"^ 
^p,d ^(p+i)d+i-p(d-i)y V 2(p-l) ) [,2nd)ir(^~^) 



p-l 2 J 



with l/r]=p{d + 2-p{d-2)). 

It is easy to relate C^^ and Kp^^. As in [16], apply (3) to fx such that 

/a(x) = /(Ax) for any x G M'^. With a := J^a |V/|2 dx, b := J^a 1/1^+^ dx, 
a := ^ + 2 — d and /3 := d^^, Inequality (3) amounts to 

Optimizing the left hand side with respect to A > shows that 

U- /'^GN\2P7 _ a + 

y-p,d ) - — — ^ ■ 

Let us consider (24). With p = 2m-i ' ^^^^ is, m = and F[u] = 
7^(p)[y] with u = /^P, it is straightforward to check that 

K- _( 2m-l \'^ d{m-mc) ^'^* _ 4 d-p{d-A) ,^1-^ J-^"'^ 

l^p.d - \, i„m ; (d+2) m.-d d(m-mi) ~ (p-1)^ <i+2-p ^^^* 

since u = always provides the equality case. Hence, using Identity (12), 
Inequality (24) amounts to 



[ \Vf\Ux+ [ \f\^+'dx-Kp,Jf in^^dx 

^ (2m-l)2 - i {m-mc) d^ (m-mc) (m-mi) {1-m)'^ 



m (1—m) * 8mKi ^ 



Using Ki = (^^_l^2')rn-d ^^*~'^ expressing everything in terms of p, we finally 



get 

Cp,d = g^(^_~^):j {{d + 2)m-d)d^ {m - nic) {m - mi) ^ — 

_ {d-p (d-4)) {d-p (rf-2)) (d+2-p (d-2)) ^ 

- 16F¥^T? a,{p) ■ 
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